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A convection-driven multiscale dynamo model is developed in the limit of low Rossby 
number for the plane layer geometry in which the gravity and rotation vectors are aligned. 
The small-scale fluctuating dynamics are described by a magnetically-modified quasi- 
geostrophic equation set, and the large-scale mean dynamics are governed by a diagnostic 
thermal wind balance. The model utilizes three timescales that respectively characterize 
the convective timescale, the large-scale magnetic evolution timescale, and the large-scale 
thermal evolution timescale. Distinct equations are derived for the cases of order one and 
low magnetic Prandtl number. It is shown that the low magnetic Prandtl number model 
is characterized by a magnetic to kinetic energy ratio that is asymptotically large, with 
ohmic dissipation dominating viscous dissipation on the large-scales. For the order one 
magnetic Prandtl number model the magnetic and kinetic energies are equipartitioned 
and both ohmic and viscous dissipation are weak on the large-scales; large-scale ohmic 
dissipation occurs in thin magnetic boundary layers adjacent to the solid boundaries. 
For both magnetic Prandtl number cases the Elsasser number is small since the Lorentz 
force does not enter the leading order force balance. The new models can be considered 
fully nonlinear, generalized versions of the dynamo model originally developed by Chil¬ 
dress and Soward [Phys. Rev. Lett., 29 , p.837, 1972]. These models may be useful for 
understanding the dynamics of convection-driven dynamos in regimes that are only just 
becoming accessible to direct numerical simulations. 


1. Introduction 

It is now generally accepted that many of the obs erved planeta r y and stellar magnetic 
fields are the result of convectively driven dynamos ( Mieschl 2005 ; Stanley fe Glatzmaier 
2010l) . Direct numerical simulations (DNS) of the complete set of governing equations 
are now routine practice, but remain limited to parameter values that are quite dis¬ 
tant from natural systems owing to the massive requirements for numerically resolving 
disparate spatiotemporal scales (e.g. Jonesl 2oTll) . The stiff character of the governing 
equations, while an impediment to DNS, provides a possible path forward for simplify¬ 
ing, or reduci ng, t he governing equa tions with the use of multiscale asymptotics (e.g. 
Julien & Kno bloch 200 7t Klein 2010t ). Balanced flows, in which two or more forces in 
the momentum equations are in balance, are particularly suitable for asymptotic anal¬ 
ysis given the subdominance of inertial accelerations. Indeed, reduced models based on 
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the geostrophic balance, in which the Coriolis and pressure gradient forces balance, have 
formed the backbone for theoretical and numerical investigatio ns o n the lar ge-scale dy¬ 
namics of the Earth’s atmosphere and oceans for over 60 years (Charnevl ll948i : [Pedlosky 


198711 . 


It is an unfortunate fact, however, that direct measurements of the forces present 
within the electrically conducting regions of natural systems are not possible. That most 
large-scale pl anetary a nd ste ll ar magne tic fields are aligned with their respective rotation 
axes (e.g. ISchubert fc Soderlundl 12011) suggests that the Coriolis force plays a key role 
in the magnetic field generation process, at least with respect to the large-scale dynam¬ 
ics. For the case of the Earth’s liquid oute r core , viscous stresses are likely to be small 
for large-scale motions there (jPozzo et al. 201 3). and observations tracking the move¬ 
ment of the geomagnetic field show th at typical con vective timescales are significantly 
longer than the rotation period (Finlay & Amit 2011). These studies suggest that the 
large-scale dynamic s within the core may be geostrophically balanced (Jault et all 1 19881 : 
Ijackson et al. 19931 iGillet et q?Jl20ld : Schaeffer fe Pais 2011 : Gillet et al. 2012 ). Alter¬ 
natively, it has been hypothesized that the Lorentz force can balance with the Coriolis 
and pressure gradient forces to y ield magnetostr ophically balanced motions at certain 
lengthscales within the core ( Roberts 1988 : Moffatt 2008h . It can be argued that only 
geostrophically balanced d ynamos have been observed with DNS (e.g. ISoderlund et al. 
2012^ King fc Buffettll2013lf (However, see the plane layer investigation of Rotvig fe Jones 
( 20021 1. 

In a seminal investigation, [Childress fc Sowardl ( 1972 ) demonstrated that laminar con¬ 
vection in a rotating plane layer geometry is capable of supporting dynamo action. They 
outlined three unique distinguished limits, or balances, in the governing equations that 
can occur based on the relative strength of the magnetic field^which they referred to as 
the weak field, intermediate field, and stro ng field limits. ISowardl ( 197411 (denoted S74 
hereafter) and Fautrelle fc Childressl ( 19821 ) (denoted FC82 hereafter) subsequently in¬ 
vestigated the weak and intermediate field limits in greater detail; for both investigations 
the flows were weakly nonlin ear. S74 s ho wed th at sta ble periodic dynamos exist in the 
weak field limit, whereas Fautrelle fc Childressl (1982) found that dynamos of interme¬ 
diate field strength may be dynamically unstable and thus lead to strong field states 
as the flow saturates. For simplicity, we refer to both the S74 and FC82 cases as the 
Chil dress-Soward dy namo model (denoted by CSDM), since they were first discussed by 
Childress fc Sowardl (1972 : and the only difference between the two models is the strength 


of the magnetic field. The basis for the CSDM is that the Ekman number Eh = v/2Q.H 2 
is taken as a small parameter, where v is the kinematic viscosity, fl is the rotation rate 
and H is the depth of the fluid layer. In this case convection is spatially anisotropic with 
large aspect ratio H U Jt> 1, wher e i characterizes the small horizontal length scale of the 
convection ( Chandrasekhai11l96ll ). By expanding the flow variables in powers of E 1 / 6 , a 
dynamical equation for the large-scale horizontal magnetic field is derived that is driven 
by the mean electromotive force (emf) generated by the small-scale fluctuating velocity 
and magnetic fields. Because the CSDM utilizes a single time scale, it applies to order 
one thermal and magnetic Prandtl numbers, Pr = v/n and Pm = v/r 7 , respectively, 
where k is the thermal diffusivity and 77 is the magnetic diffusivity of the fluid. 

The CSDM continues to provide an invaluable tool for understanding the differences be¬ 
tween large-scale and small-scale dynamo action, as well as aiding the interpretation of re¬ 
sults obtained from DNS studies. S74 showed that a dynamo driven by small-scale rapidly 
rotating convection is necessarily large-scale. Many of the key features predicted by S74, 
such as the vertical structure and oscillatory temporal evolution of the large-scale mag- 
netic field, have been confirmed with low Ekman number DNS (e.g. St ellm ach fe Hansen 
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20041 ). Recent work has employed t he CSDM for investiga ti ng kinematic dynamo ac- 


tion driven by rotating convection ( Favier fc Proctorl 1 201 .'ll) . iMizerski fc Tobiah ( 20131 ) 
have extended the CSDM to include the effects of fluid compressibility via the anelastic 
approximation, showing that density stratification tends to delay the onset of dynamo 
action and reduces the strength of the resulting magnetic field. 

In the present work we develop a new multiscale dynamo model that possesses many 
similarities with the CSDM, but also significant differences. We refer to the new model 
as the quasi-geostrophic dynamo model, or QGDM. Like the CSDM, the small-scale, 
convective dynamics are geostrophically balanced to leading order, but remain time- 
dependent and fully nonlinear in the QGDM; our small-scale dynamical equations are 
a magnetically modifi ed version of the quasi-geostr ophic convection equations developed 
bv ljulien et all 1 19981) (see also Julien et al. ( 20061 ) for a generalized development). The 
strength of the magnetic field in the QGDM is asymptotically larger than that considered 
by FC82, but is not a strong field in the sense that the Lorentz force does not enter the 
leading order force balance. In addition to the small horizontal convective lengthscale t, 
the QGDM includes large-scale horizontal modulations Lx that allow for a diagnostic 
thermal wind balance in the large-scale momentum equations. Coupled with the large- 
scale heat equation, the large-scale model (in the absence of a magnetic field) is the 
equivalent of the planetary geostrophic equa t ions commonly emp loyed in oceanography 
( Robinson fe Stommel 1959 ^felajidgi[ 1959: iGrooms et a/.ll201l|) and atmospheric sci¬ 
ence ( Phillipdll963 : Dolaptchiev fc Kleinll2009l) . (The so-called /3-effect has been neglected 
in the present work though it can easily be incorporated into the current model, e.g. see 


Grooms et al\ ( 201ll) ). Because of the additional scale Lx, the QGDM also possesses a 


non-zero large-scale vertical magnetic field that is not present in the CSDM. Moreover, 
the new model utilizes three disparate timescales characterizing the small-scale convective 
dynamics, the large-scale magnetic evolution timescale, and the large-scale thermal evo¬ 
lution timescale. Both low and order one magnetic Prandtl number cases are considered. 
The QGDM can be simulated numerically to gain insight into large-scale magnetic field 
generation in planets and stars , where we can utilize t he success of previous non-magnetic 
work (e.g. Sprague et al. 200fil : Julien et~. I 2012ol f&l: iRubio et a/.l 120 1 ll : Stellmach et al 


2014!), and the increases in computational power that have occurred since the work of 


CS72. 

In section [2] we present the asymptotic development of the model. In section [3] we 
discuss some of the important features of the QGDM, its relationship with the CSDM, 
and possible applications to natural dynamos. Concluding remarks are given in section 

SI 


2. Model Development 

2.1. Governing Equations 

We consider a rotating plane layer geometry with rotation vector Cl = Liz and constant 
gravity vector g = — gz. The horizontal boundaries are located a vertical distance H 
apart with the layer heated from the bottom boundary and cooled at the top boundary. 
For simplicity, the fluid is assumed to be Newtonian and Boussinesq with density p and 
thermal expansion coefficient a. The governing equations are then non-dimensionalized 
utilizing the scales U , i, IjlA , AT, V and B for the velocity, length, time, temperature, 
pressure, and magnetic field, respectively. In the rotating reference frame with coordinates 
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( x , y , z) the dimensionless governing equations are then 


d t u + u • Vu H- z x u = —Eu S7p + M B • VB + T 8 z H-V 2 u, (2.1) 

Ro Re 

d t 0 + V-(u 0) = -j-V 2 0, (2.2) 

Pe 

9 f B = Vx(uxB) + tt-V 2 B, (2.3) 

Rm 

Vu = 0, VB = 0, (2.4) 


where the velocity, pressure, temperature, and magnetic field are denoted by u = (it, v, w), 
p , 9 , and B = (B( x \ B^ y \ B^^j, respectively. Both the hydrostatic centrifugal force and 
magnetic pressure have been absorbed in the pressure gradient term Vp. The dimension¬ 
less parameters are defined by 


u 

r U 2 ’ R 

Ui 

5 

V 

Ut 

Pe = —, 

K 

V 

Eu= PIP’ 

(2.5) 


M B2 

Rm 

Ui 


(2.6) 


ppW 1 

V ' 



Here, Ro is the Rossby number, T is the buoyancy number, Re is the Reynolds number, 
Pe is the Peclet number, Eu is the Euler number, M is the ratio of magnetic energy to 
kinetic energy, and Rm is the magnetic Reynolds number. The permeability of free space 
is denoted by p. We recall that the thermal and magnetic Prandtl numbers characterizing 
the fluid properties are related to the above parameters via the relationships Pe = RePr 
and Rm = RePm. Importantly, in the present work we denote dimensionless parameters 
that are based on the large-scale depth of the fluid layer with a subscript H, whereas 
parameters without this distinction utilize the small convective scale t. 

In what follows we assume that the horizontal bounding surfaces are impenetrable, 
stress-free and perfect thermal and electrical conductors such that 

d z u = d z v = w = d z B = d z BW = B W = 0, at z = 0, H/i, (2.7) 


9=1 at z = 0, and 0 = 0 at z = H/i. (2-8) 

The development given here parallels that of Julien fe Knoblochl ( 2007 ) and Grooms et al\ 
<Iml) for other multiscale models. We find that three disparate timescales are neces¬ 
sary to completely describe the temporal evolution of the dynamo; these are (1) the 
timescale of the small-scale convective fluctuations (denoted by t), (2) the large-scale 
magnetic evolution timescale (r), and (3) the large-scale thermal evolution timescale 
(T). The small-scale, rapidly varying spatial and temporal coordinates are denoted as 
(x, t), and the slowly varying large-scale coordinates are given by (X, r, T). The fields 
are decomposed into mean and fluctuating components according to 


u(x, X, t, r, T) = u(X, r, T) + u'(x, X, t, r, T), (2.9) 

p(x, X, t, r, T) = p(X, r, T) + p'(x, X, t, r, T), (2.10) 

0(x, X, t, r, T) = 0(X, r, T) + 0'(x, X, f, r, T), (2.11) 

B(x, X, t, t, T) = B(X, r, T) + b'(x, X, t, r, T), (2.12) 


and the fluctuating magnetic field vector is b' = ( b'^ x \ b' <y \ b'^). In Figure[l]a schematic 
illustrating the general features of the model is shown. The mean quantities are defined 
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Figure 1: Schematic summarizing the general features of the multiscale dynamo model. 
The large-scale coordinate system is X = (X , Y, Z) and the small-scale coordinate system 
is x = ( x,y,z ); the shaded boxes depict small-scale subdomains embedded within the 
single large-scale domain. Large-scale mean dependent variables are given by (u, p, 6 , B) 
and the small-scale fluctuating dependent variables are (u 1 ,p' ,9' ,h'). Distinct equation 
sets are developed for the large-scale and small-scale domains, with the mean dynam¬ 
ics controlling the background state for the fluctuating dynamics and the fluctuating 
dynamics feeding back onto the mean state via small-scale fluxes. 


as averages taken over the fast scales according to 




't',V 


/(x,X,t,r,T)dx£ft, /' = o, 


(2.13) 


with the small-scale fluid volume denoted by V. With multiple scales the differentials 
become 


d t ->d t + -j-e? T + -7-dr, 

J±r At 

v -> v + v, 


(2.14) 

(2.15) 


where now V only acts on the small-scale coordinates x and the large-scale, mean gradient 
operator is defined by 


V = (—ft,- —r)v —I 


\A x dx '7^ dv — 2 dz ^ 


The various aspect ratios are defined according to 


. . Lx H t t 

Ax = Ay = -, Az — —, A r = —, Ay = —. 

i z i t 1 T 


(2.16) 


(2.17) 


For simplicity we assume isotropic horizontal large-scale modulations (i.e. Ax = Ay) in 
the present work, though the approach remains completely general. 

The governing equations become 


Z X u = 


(dt + J-d T + U + V- K) + V- K) + 

— Eu (V + V) p + M (B • V + B • V) B + T 9z + (v 2 + V 2 ) u, 


1 

Ro 


d t + -y-d T + ——dy ) 9 + V- (0u) + V- (0u) = 
At Ay 


Re 

1 

T'e 


(v 2 + v 2 ) 


(2.18) 


(2.19) 
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(dt + A-dr + B = V x (u x B) + V x (u x B) + -^ fv 2 + V 2 ) B, (2.20) 

\ A r At J ixttl V / 

(V + V) • u = 0, (V + V) • B = 0. (2.21) 

The superscript “tr” appearing on some of the velocity vectors denotes a transpose. 

Averaging the equations over fast temporal and spatial scales (x, t) results in the mean 
equations 


(A-d T + A-dr) u + V- (u(u) tr ) + V-(u'(u') tr ) + A- z x u = 
\ A r At J Ro 

— — — ^ 1 — 2 _ 

— Eu\7p + MFl + TOz + —— V u. 

Re 

(r dr + j^ dr ) ° ++ = 7 ^ 0 , 

A-d r + —9t ) B = Vx(uxB) + Vx (u' x b') + A- V 2 B, 
,A t A t J Rm 


( 2 . 22 ) 


(2.23) 

(2.24) 


V ■ u = 0, V • B = 0, (2.25) 

where the mean Lorentz force isFi = B-VB + b' • Vb'. 

The fluctuating equations are found by subtracting the mean equations from the full 
equations to give 

(dt + A-dr + A-dr) u' + u-Vu' + u'-Vu-V- (u'(u') tr ) + A- z x u' = 

V A r A t ) Ro (2 . 26) 

-Eu (V + V)p' + MF' L +r0'z + A- (V 2 + V 2 ) u', 

(d t + A-dr + A~d T ) & + u' • V 0 + u • V6»' + V • (u'0') - V • (VT) = 

V At a t ) (2.27) 

T( v » + v 2 )»', 

( a * + + i*) b ' = 

Vx(uxb') + Vx(u'xB) + Vx(u'xb')-Vx (u' x b')+ ( 2 - 28 ) 

V x (u x b') + V x (u' x B) + V x (u' x b') + — fv 2 + V 2 ") b', 

Rm \ J 

(V + V) • u' = 0, (V + V) • b' = 0, (2.29) 

where the fluctuating Lorentz force is = B • Vb' + b' • V B + b' • Vb' — b' - Vb'. 


2.2. Asymptotics 

We note that up to this point, no approximations have been made; the equations have 
simply been split up into mean and fluctuating components. In the present work we are 
interested in the development of a multiscale dynamo model that preserves geostrophic 
balance on the small (fluctuating) scales. Because the fluctuating and mean dynamics are 
coupled, a posteriori we find that the large-scales are required to also be geostrophically 
balanced in the horizontal directions, and hydrostatic in the vertical direction, resulting in 
a thermal wind balance. The relevant asymptotic limit for the present model is therefore 
the small-Rossby number limit, Ro = e —>• 0. 
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When constructing an asymptotic model with more than one small (or large) parame¬ 
ter, it is necessary to establish the so-called distinguished limits, or the asymptotic scaling 
relationships between the various parameters. In the present context this implies that we 
must relate, in an order of magnitude sense, both the non-dimensional parameters intro¬ 
duced in equations (12.5I) - (12.6I) . as well as the aspect ratios defined by equations (12.17|) . to 
the Rossby number e. This procedure is well known, and hinges on the identification of 
dominant balances in the governing equations that will allow for a mathematically and 
physically meanin gful reduced model; by definition the procedure is circular in nature 
(e.g. see lBender fe Orszad 2010h . In this regard, t o ensure geostr o phica lly balanced con¬ 
vection on the small scales, we follow section 2 of Sprague et al\ ( 2006h and employ the 
following distinguished limits 


A? = 


At 


Eu • 2 , r = e” 1 T, Re = 0(1), Pe = 0(1), (2.30) 


where the reduced buoyancy number T = 0(1). We note that with these scalings, the 
small-scale velocity field u' is 0(1). We further assume the large-scale pressure p , mag¬ 
netic field B and temperature 9 to also be 0(1) in magnitude. 

The dimensional scales employed in the dimensionless parameters given by equations 
(ESD-lEH) will be those of the largest amplitude. For convenience the parameters Eu, T, 
and M will therefore be based on ‘mixed’ scales in our asymptotic development in the 
sense that the pressure, temperature, and magnetic field amplitudes are based on the 
large-scale quantities (p, B, 9) whereas the velocity scale is based on u'. In this sense, the 
dimensional magnitudes of these variables is absorbed in the various non-dimensional 
parameters. In section 13.41 we discuss the renormalization of these non-dimensional pa¬ 
rameters such that they are based on equivalent (large or small) dimensional scales. 

It now remains to determine the distinguished limits of Ax, A r , M , and Rm. To pro¬ 
vide some explanation as to how these limits are achieved we examine various terms in 
the governing equations that we wish to retain in the final reduced model. The distin¬ 
guished limit of Ax can be determined by restricting the large-scale horizontal motions 
to be geostrophically balanced; utilizing equation (12.221) we have 

-ZXU« — Vxp, (2-31) 

e e-A x 


where Vx = (dx,dy, 0). We then get 


u ~ (eA x ) 1 ■ 


(2.32) 


Turning our attention to the mean temperature equation (12.23[) . we wish to keep hori¬ 
zontal advection of the mean temperature by the mean velocity field at the same order 
as the temporal evolution of the mean temperature so that 

e 2 d T 9 ~ u • Va -9, (2.33) 

Ax 

which leads to 

u ~ e 2 A x . (2.34) 

Combining (12.321) with (12.341) then shows that we must have 

Ax = e~ 3/2 , (2.35) 

and so u = 0(e 1 ^ 2 ). We note that this also implies large-scale horizontal gradients 
are smaller than the corresponding large-scale vertical gradients by a factor Az/Ax = 

0(A/ 2 ). 
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The distinguished limits of A T , M, and Rm are determined by examining the influence 
of the magnetic field. To generate dynamo action on the large scales via coupling with 
the fluctuating dynamics, we necessarily require the presence of the second term on the 
righthand side of (12.241). where (u' x b') is referred t o as the mean electromotive force, 
or me an emf (e.g. |Parkeijll955l : ISteenbeck et al. 1966t Steenbeck fe Krausell 19661 iMoffatt 
Il978l) . Utilizing the above scalings for Ax and u, and retaining only the largest terms, 
the order of magnitude of each term present in the mean induction equation is then 


e 3 / 2 B : eb' 


-B, 


(2.36) 


B 

A t Rm 

w here we recall that the fluctuating velocity field is order one, consistent with lSpraeue et all 
(2006]). The colons used above are meant to signify that we are making a comparison 
of the magnitude of each term in the original governing equation and we have retained 
B = 0(1) for clarity in identifying the various terms. To retain time dependence of the 
mean magnetic field, we require that the first and third terms be of the same magnitude 
such that 

B ~ eA T b'. (2.37) 

By noting that all the large-scale gradients are small relative to the small-scale gradients, 
the order of magnitude of the four largest terms (the first on the lefthand side, and the 
fifth, sixth and eighth terms on the righthand side) in the fluctuating induction equation 
(|2.28l) are 

b' : B : b' : (2.38) 

Rm 

We require that ohmic dissipation is present in the final model, at least with respect 
to the small scales. The largest term present above would then be that associated with 
stretching the mean magnetic field such that 


B 


b' 

Rm 1 


(2.39) 


i.e. b' ~ Rm. Finally, because we’re investigating dynamo action , we r equire that_the 


Lorentz force enter the small scale momentum dynamics. From ISprague et al. 
and examination of equation (12.26[) , we know that this requires 

MBb' = 0(1). 

Combining this scaling with (12.391) yields 

M = 0(Rm~ 1 ). 


( 120061 ) 


(2.40) 


(2.41) 


We can now determine the size of Rm by returning to the mean induction equation. For 
large-scale dynamo action, magnetic diffusion cannot dominate the mean emf so that 
comparison of the third and fourth terms, respectively, given in (12.361) yields 

e 2 _ 

eb' > -B. (2.42) 

Rm 

Using the above scaling with (12.391) shows that 

Rm > e 1/2 , (2.43) 

with the lower bound being the case Rm = 0(e 1 / 2 ). The particular distinguished limit 
taken for Rm will in turn determine the size of the magnetic Prandtl number since 
Rm = RePm. Taking Rm = 0(e 1 ^ 2 ) corresponds to the low Pm limit since Re = 0(1). 
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An alternative scaling would be to take Rm = 0(1) and thus Pm = 0(1). The main 
focus of the present manuscript is for the Rm = C^e 1 / 2 ) limit given that it ties directly to 
the CSDM and is implicitly low Pm. We also present the form of the governing equations 
for the Rm = 0(1) case, but note additional complications that arise in this limit which 
will require future analysis that is beyond the scope of the present work. 

2.2.1. The Rm = O^ 1 / 2 ) Limit 

Here we take Rm = e^^Rm, where we define the reduced magnetic Reynolds num¬ 
ber Rm = 0(1). Returning to relationships (12.371) and (12.411) we can now define the 
distinguished limits of A r and M as 

A t = e“ 3/2 , M = e“ 1/2 , (2.44) 

showing that three disparate timescales are now present in the QGDM. 

With our distinguished limits defined, we now follow CS72 and expand all variables in 
powers of e 1 / 2 , e.g. 

u = Uo + Ug + e 1 /" + u i/ 2 ^ + e ( u i + u i) + • • • ■ (2.45) 

Plugging the perturbation expansions into the governing equations and collecting terms 
of equal magnitude, the leading order solenoidal conditions on the mean velocity and 
magnetic fields gives 

d z u J 0 = 0, d z Bf = 0. (2.46) 

Owing to our use of impenetrable, perfectly conducting boundary conditions we require 
w o = B o = 0. At the next order the solenoidal conditions yield 

V.v-uo + dzvJi/2 = 0, Vx-Bo + dzB^/ 2 = 0. (2.47) 

At order 0(e _1 ) the horizontal mean momentum equation gives 

z x u 0 = 0, (2.48) 

showing that uo = 0. Utilizing continuity it follows that uJi / 2 = 0. The next three orders 
of the mean horizontal momentum equation are then geostrophically balanced, 

1 3 

z x Ui = -VxPi— 1 /2 > i= 2 ,lj 2' < ' 2A9 ’ > 

Carrying the expansion out to O(e) shows that u 2 is magnetostrophically balanced, 

z x u 2 = — VXP 3/2 + B 0 • VxB 0 + B[/ 2 dzBo. (2.50) 

The vertical mean momentum equation yields hydrostatic balance for the first four orders, 

dzPi= f0i, z = 0, |, 1, |. (2.51) 

Taking the curl of equation (12.49[) shows that 

Vx ■ Ui = 0 , i = 0 ,^, 1 ,|, (2.52) 

and thus uJ* = 0 for i = 0,... 2. Combining equation (12.491) for i = 1/2 and (12.511) for 
i= 0 we then obtain the well-known thermal wind relations 

dzv 1/2 = TdxO 0 , dzui/2 = —Td y 9 q. 


(2.53) 
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Because the large-scale velocity field is horizontally divergence free we can define the 
large-scale geostrophic stream function as '& 1/2 =Po such that 


U1/2 = -V X 'f'i/ 2 Z. 


(2.54) 


We mention that the mathematical structure of the thermal wind relations (12.5311 yields 
any barotropic (i.e. depth invariant) geostrophic large-scale flow undetermined given that 
U 1/2 only appears with first-order Z-derivatives. In some instances this barotropic flow 
plays a signific ant dynamical role a nd must be consistently determined via an evolution 
equation (c.f. iDolaptchiev fc Klein 2009 ). However, in the present work the evolution 
equation for such a mode enters at a much higher, subdominant order, and we do not 
consider it any further. 

Proceeding to the mean horizontal induction equation, at O(e) we have 


0 = dz 


zx (u' 0 x bg) 


(2.55) 


which, in general, requires that bg = 0. Alternatively, from the scaling relationship given 
by (12.391) . we require that b' = 0(e 1 / 2 ) if B = 0(1). At 0(e 3 / 2 ) we get 


d r B 0 = d z 


z x 


( U 0 X b i/2 


Km 


(2.56) 


l _ (.%') _ (y) 

where B 0 = (B Q , B 0 ,0). The leading order vertical component of the mean induction 
equation becomes 


8 t b[% = d x 


y • 


( U 0 X b i /2 


~dy 


( u o x b i/ 2 ) 


1 a2 -nC) 


Rm 


dzB\y 2 , (2.57) 


where we see that a non-trivial mean vertical magnetic field requires the presence of a 
large-scale horizontal modulation in the sense that B{/ 2 = 0 if dx = dy = 0 (this fact 
can also be seen from equation (12.4711 1. The fluctuating induction equation then gives 


0 = Bn • Vj.uA + ^V 2 b' 


Rm 

where Vj. = {d x , d y , 0). 

The leading order mean temperature equation gives 


1/25 


dz(w' 0 6' 0 ) = 0, 

showing that 9' 0 = 0. At the next order we have 

d T do + dz 


'( w o d 'i/2) 


= 0, 


again showing that d' x , 2 = 0 and thus d T 0 0 = 0. Finally, we have 


drO 1/2 + drOo + u ±/ 2 ■ S 7 x @0 + dz(vjQ 0 i) — ~p^dz® o- 


(2.58) 

(2.59) 

(2.60) 

(2.61) 


To avoid secular growth of the mean temperature on the timescale r we set <9 r $i/2 = 0; 
the mean heat equation then becomes 


drO 0 + u 1( / 2 • + dziw^O'Y) = j^d z 9o. 


(2.62) 


Alternatively, one can average equation (12.611) over the timescale r to obtain an equation 
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identical to (12.6211 . with the exception that the averages must then be interpreted as 
occuring over (x, t, r). At O(e) the fluctuating temperature equation is 

1 


d t e[ 


* 0 -v±o[ + w , 0 d z 6 0 = 1 ^v 2 e' 1 . 


(2.63) 


At 0(e 2 ) and 0(e 3 / 2 ) the fluctuating momentum equation yields, respectively 


Vp'=0, i = 0, —, 


showing that p' 0 = p[, 2 = 0. Geostrophy occurs at 0(e 3 ) and 0(e 1 / 2 ) 


z x u' = -Vp' 


2+1 ’ 


i = 0 


Additionally, mass conservation at 0(1) and C^e 1 / 2 ) gives 


V-u' = 0, 


* = °’ 2 ’ 


(2.64) 


(2.65) 


( 2 . 66 ) 


which, along with equations (12.651) . yields the Proudman-Taylor theorem acting over the 
small vertical scale z 

dz (u',p' +1 ) =0, i = 0,i. (2.67) 

It is natural to wonder why the Proudman-Taylor theorem is satisfied over the small 
scale z and not the axial domain scale Z. It is well known from the linear theory of 
rapidly rotating convection that the leading order, geostrophically balanced convection 
possesses order one variations over the height of the fluid layer (our large-scale coordi¬ 
nate Z) ( Chandrasekhaifl961 ). The only way to be consistent with the Proudman-Taylor 
theorem, whilst allowing convection to occur on the axial domain scale, is for it to be 
satisfied over the small vertical scale z. We recall that this is the same scale as the horizon¬ 


tal convective scales {x,y). This fact appears to be first noted by IStewartson & Cheng 
( 19791) : apart from the presence of viscous and body forces in the present work, their 
vorticity (2.7b) and vertical momentum (2.7a) equations are identical to our small-scale 
counterparts given by (12.6911 and (12.701) . 

The prognostic momentum equation then appears at 0(1) 


<9 4 Uq + Uq • Vj_Ug + z x u' x = -Vp 2 - d z p’i z + TO^z, + B 0 • V+b) 


1/2 




( 2 . 68 ) 


where we note the key distinction with the CSDM is the simultaneous appearance of 
both the advection and Lorentz force terms. In light of equation (12.651) . we can define the 
small-scale geostrophic stream function as i/j' 0 = p[ such that the fluctuat ing horizontal 
velocity field is given by (u' 0 ,v' 0 ) = (—dyil’b,d x ip' 0 ). As noted bv ISpraeue et al\ ( 20061 ). 
equation (12.681) still depends upon the small-scale vertical scale z; to avoid secular growth 
on this scale it is necessary to impose solvability conditions. These solvability conditions 
amount to operating on equation (12.681) with z • Vx and z-, respecti vely, and averaging 
over the small vertical scale z to obtain (see also lCalkins et aill2013l) 


a t VlV'o + ^o.Vi^)-az«/ 0 =2-Vx(Bo-Vx(b' 1/2 )) + ^Vi^, (2.69) 

d t w' 0 + JM,w , 0 ) + dz% = f(e' 1 )+z. (Bo-V+<b' 1/2 ))+-i-Vi^, (2.70) 

where J(F, G) = d x Fd y G — d x Gd y F , the angled brackets denote a spatial average over 
z, and the vertical vorticity is (o = ViV’o- 
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Upon rescaling the velocity with the small-scale viscous diffusion time such that U = 
v/L, the closed set of reduced equations is given by 


U1/2 = -V x ^i/ 2 z, 

_ Ra-rr 

dz'S’i /2 - -p^O o, 


1 q2 r 


d T do + ui/2 • V x 0o + dz{w' 0 6' 1 ) = —d 2 z 6 0 , 


9 r B 0 = d z 


z x 


( u o x < b 'l/ 2 >) 


Fm 


o t b[% = dx 


9 ■ ( u o x < b i /2 >) -°y x • ( u o x < b i/ 2 >) 


1 d\B { 


Pm 


i/ 2 ) 


d x B ( 0 X) 


dvB ( 0 v) + d z B[% = 0, 


S t ViV>0 + J(i/> 0, ViV'o) - dzw'o = z • V x (b 0 . V x (b' 1/2 )) + 


dtw o + J(ip o, w' 0 ) + dz^o = ^(0[) 


Bn 


vM/1) + vv 0 , 


W + J(^o, (0i)) + 

0 = B^-V ± u' iivKbi/j), 

Fm 

dx( b i/2) + d y( b ijl) = °- 


(2.71) 

(2.72) 

(2.73) 

(2.74) 

(2.75) 

(2.76) 

(2.77) 

(2.78) 

(2.79) 

(2.80) 
(2.81) 


The reduced Rayleigh number, consistent with the linear theory of rapidly rotating 
convection ( Chandrasekhar-! fl961 ). is defined by Ra = e^Ran = E^ 3 Ran, where the 
Rayleigh number is given by 

gaATH 3 


Roh = 


VK 


(2.82) 


In addition,^the reduced magnetic Prandtl number appearing in the above system is 
denoted by Pm. 

The final reduced system is 10th order in the large-scale vertical coordinate Z. Written 
in terms of the reduced variables, the ten boundary conditions become 


do = 1 at Z = 0, do = 0 at Z = 1, 


d z P>o = b i /2 =w' o =0 at Z = 0,1. (2.84) 

Although the details of obtaining numerical solutions to the reduced system are beyond 
the scope of the present work, we mention that several approaches have been developed to 
handle both multiple scales in space and time. In particular, so-called heterogene ous m ul¬ 


(2.83) 


tiscal e methods (HMM) have been successfully applied to numerous applications (IE e t al 


200711 . For example, a modifed form of HMM has recently been de velo ped for simulating a 


multi scale model of upper ocean Lan g muir circulation ( Chini et al. 20091: Malecha et al. I 
203). In addition, Haut fc Wingatel ( 2014h have developed a time-stepping algorithm 


for problems with temporal scale separation which may be of use for the present model. 


2.2.2. The Rm = 0(1) Limit 

The asymptotic development of the Rm = 0(1) parallels that for the Rm = 0(e 1 ^ 2 ) 
case given above. For this reason we omit many of the details and focus only on the 
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differences between the two limiting cases. We emphasize again that the Rm = 0(1) 
limit corresponds to Pm = 0(1) since Re = 0(1). Consideration of relationships (12.371) 
and (12.411) with Rm = 0(1) shows that we must have 

A T =e~\ M = 0(1). (2.85) 


In addition, relationship (12.391) shows that the mean and fluctuating magnetic fields are 
now of the same order, i.e. B ~ b'. The main difference that occurs in the final reduced 
model for the Rm = 0(1) case is the form of the induction equations, which at leading 
order become 


d T B Q =d z 


j(*) 


d T B 1/2 = d x y • (uq x (bo)) 


z x (uo x (b^)) 

-dy 


K X (bg)) 


a t (b') + u' • Vj_(bo) = B 0 • V_lUq + —Vi(b'). 


( 2 . 86 ) 

(2.87) 

( 2 . 88 ) 


The fluctuating equations then take a mathematically equivalent form of equations (I2.77D - 
(12.79[) and (12.811) by making the substitution (b^ 2 ) —> (b(,). 

The boundary conditions for the Rm = 0(1) case are identical to those given by (12.831) - 
(I2.84D . Equations (I2.86D - Q2.87D shows that due to the absence of ^-derivatives, no mean 
magnetic field boundary conditions can be satisified without the inclusion of magnetic 
boundary layers. Future work is necessary to examine the effects of these boundary layers. 


3. Discussion 

3.1. The Rm = (^(e 1 / 2 ) limit 

Here we reiterate some of the key features of the low magnetic Prandtl number QGDM 
and interpret the various terms in the equation set (12.71D - (12.811) . Equations (I2.71l) - (12.72l) 
are statements of geostrophic balance and hydrostatic balance on the large horizontal 
and vertical scales, respectively. Importantly, these equations are diagnostic since they 
contain no information about the temporal evolution; this is a well-known characteristic 
and the prognostic dynamics are obtained from the mean heat equation (12.731) . From 
equation (12.711) we see that viscous diffusion is negligible on the large scales. The mean 
heat equation given by (12.731) shows that the mean velocity field is strong enough to allow 
advection of the mean temperature over the large horizontal scales (X, Y) as shown by 
the second term on the lefthand side, whereas the third term represents the influence of 
convective feedback from the fluctuating scales back onto the large-scales. Additionally, 
the presence of large-scale thermal diffusion over the vertical dimension is represented 
by the term on the righthand side of equation (12.731) . 

Both of the mean induction equations given by (12.741) and (12.75[) contain time depen¬ 
dence, the mean emf that represents the feedback from fluctuating velocity and magnetic 
field dynamics, and ohmic diffusion on the large vertical scale Z. Thus, ohmic dissipation 
dominates viscous dissipation on the large scales since, at this order of approximation, 
viscous diffusion is not present in equations (I2.71I) - (I2.72I) . This feature is consistent with 
the Pm <C 1 limit considered here. Moreover, the QGDM is characterized by a horizontal 
magnetic field that is 0(i?m _1 ) stronger than the vertical magnetic field. 

Equations (12.77[) - (12.79l) are a magnetic versi o n of the quasi-geostrophic co nvection 
equations originally developed bv ISpraeue et al. ( 20061) and Julien et al. ( 20061 ). An im¬ 
portant feature of equations (12.771) - (12.781) is that both advection and diffusion only occur 
over the small horizontal scales ( x , y ) due to the anisotropic spatial structure of low 
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Rossby number convection. The fluctuating vorticity equation (12.771) contains time de¬ 
pendence, advection on the small-scales, vortex stretching is represented by the third 
term on the lefthand side, with the remaining terms on the righthand side being the 
Lorentz force and viscous diffusion, respectively. The fluctuating momentum equation 
is given by (12.781) and, like the fluctuating vorticity equation, contains time dependence 
and horizontal advection. The vertical pressure gradient is given by the third term on 
the lefthand side of equation (12.781) . with the fluctuating geostrophic stream function 
tpQ acting as pressure since the flow is geostrophically balanced; the remaining terms on 
the righthand side are the buoyancy force, the Lorentz force, and viscous diffusion. The 
fluctuating heat equation (12.791) is also characterized by time dependence and horizontal 
advection and diffusion, with the third term on the lefthand side representing advection 
of the mean heat by the fluctuating vertical velocity over the large-scale Z. 


The fluctuating induction equation (12.801) is identical to the equation given in S74, 
where he noted that the fluctuating magnetic field is induced by stretching the mean 
magnetic field with the small-scale strain Vj_Uq (consistent with our low Rm approxima¬ 
tion). The absence of the time derivative shows that the fluctuating magnetic field adjusts 
instantaneously relative to the fluctuating convection, showing that magnetic diffusion is 
much more rapid than momentum diffusion on the small-scales; this effect is consistent 
with a small magnetic Prandtl number limit. Alfven waves are therefore damped out on 
the small horizontal scales. As with the fluctuating vorticity, momentum, and heat equa¬ 
tions, only horizontal diffusion is present in the fluctuating induction owing to spatial 
anisotropy. Taken with equations (|2.77jl - (|2.79j) . we see that both viscous dissipation and 
ohmic dissipation are important features for the small-scale dynamics. The small-scale 
induction equation also shows that both the present model and the CSDM require the 
presence of a non-trivial mean magnetic field for the development of a dynamo; this 
shows that the resulting dynamo is therefore large-scale and typical of low Rm dynamos 
dlVToffattlh 978li . 


3.2. The Rm = 0(1) limit 


The Rm = 0(1) case is associated with an order one magnetic Prandtl number. For 
this reason, the Pm = 0(1) QGDM may be particularly important for relating to DNS 
studies where reducing Pm to physically realistic values is impossible given the mod¬ 
est Reynolds numbers attainable with current computational resources. Given that the 
energy and momentum equations are identical to the low Pm QGDM discussed in the 
previous section, we focus here on the differences associated with the form of the induc¬ 
tion equations (12.861) - (12.881) . 


The distinguished limit A T = e -1 shows that, in comparison to the Rm = 0(e 1//2 ) case, 
the mean magnetic held now evolves on a faster timescale. The M = 0(1) limit indicates 
that the magnetic and kinetic energies are now equipartitio ned w hen Pm = 0(1) — a 
result that appears to be consistent with DNS investigations (jStellmach fe Hansen 200 fj. 

The absence of ^-derivatives in equations (12.86l) - (12.87l) shows that ohmic dissipation 
is weak on the large-scales and limited to magnetic boundary layers adjacent to the 
horizontal boundaries. The small-scale induction equation (12.881) contains time depen¬ 
dence, advection by the fluctuating velocity held, along with stretching and diffusion. 
The presence of (9 t (t>Q) results in a fluctuating magnetic held that no longer adjusts in¬ 
stantaneously to the fluctuating velocity held as it does when Pm <C 1. For this case, 
Alfven waves are present on the small fluctuating scales. 
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Parameter | 

Present Work | 

Soward | 

Fautrelle & Childress 

Eu 

e" 2 

e" 2 

e" 2 

Re 

0(1) 

e 1/2 

e V 2 

M 

e -1/2 , 0(1) 

£ 3/2 

e 1 / 2 

r 

Pe" 1 

Pe -1 

Pe" 1 

Pe 

0(1) 

e 1 / 2 

e 1/2 

Rm 

e 1/2 , 0(1) 

e 1/2 

e 1 / 2 

Pr 

0(1) 

0(1) 

0(1) 

Pm 

e 1/2 ,0(l) 

0(1) 

0(1) 

Ax 

£ —3/2 

- 

- 

Az 

e" 1 

e" 1 

e" 1 

a t 

£ —3/2 

£— 3/2 

e -3/2 

Ax 

e" 2 

— 

— 


Table 1: Co mparison betw een t he different distinguished li mits taken in the present work 
and those of Soward ([19741 ) and Fautrelle & Childress 198211 . relative to equations (12.221) - 
(12.29f) . We note that both Soward and Fautrelle & Childress considered weakly nonlinear 
motions and only a single (slow) timescale and a single large-scale coordinate in the 
governing equations such that d t = (?t = dx = &y = 0- In a ll models the small parameter 
is the Rossby number, Ro = e, and the reduced buoyancy number F = 0(1). 


3.3. Relationship with the Childress-Soward dynamo model (CSDM) 

In Table [I] we show the different distinguished limits taken between the present work and 
those taken by S74 and FC82. Although we have considered two limits of Rm that result 
in different limits of M, Pm, and A T , the present discussion will be focused on the Rm = 
0(e 1 / 2 ) (low Pm) QGDM. Both the weak and intermediate field forms of the CSDM can 
be derived directly from equations (I2.22I) - (I2.29I) by employing the distinguished limits 
listed under the columns labeled “Soward” and “Fautrelle & Childress”, respectively, 
and taking d t = dr = dx = dy = 0. The consequence of taking Re = Pe = 0(e 1 / 2 ) is 
that the resulting models are weakly nonlinear since (horizontal) viscous diffusion enters 
at higher order than the advective nonlinearities in the fluctuating momentum equation. 
While in the present work we have considered both low and order one Pm, the CSDM 
employs Pm = 0(1). Additionally, the magnetic field scaling M s = e 3 / 2 employed by 
S74 and M FC = e 1 / 2 employed by FC82 contrasts sharply with our M = e -1 / 2 limit; 
the result is that magnetic energy is significantly larger than the kinetic energy in the 
present work. 

CS72 described three different classes of dynamos, distinguished by the relative strength 
of the magnetic field. In their work, the classification was based on the magnitude of the 
large-scale Hartmann number defined as 


Ha H 


BH 

C hpvrtf 1 2 


(3.1) 


In the CS72 terminology, the weak field, intermediate field, and strong Held dynamos are 
characterized by Hartmann numbers of magnitude Ha^ = 0(1), Ha l H = 0(e _1,/2 ) and 
Ha s H e -1 / 2 , respectively. FC82 more precisely identified the strong field regime by the 
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scaling Ha s H = e -3 / 2 since this results in a Lorentz force that is comparable in strength 
to the Coriolis force. By employing our distinguished limits, the large-scale Hartmann 
number in the present work is given by 


Ha* H = 


(MRm) 1/2 


(3.2) 


Thus, by the classifications of CS72 and FC82 the present model considers magnetic 
fields with strengths that lie between the intermediate and strong field limits. As the 
fluctuating momentum equation (12.681) shows, both the Lorentz force and horizontal 
viscous diffusion enter the prognostic equation at the same order in our model; the result 
is that the small-scale Hartmann number is unity, Ha* = 1. 

A strong field dynamo is synonomous with the magnetostrophic balance, whereby the 
Coriolis, pressure gradient and Lorentz forces are of comparable magnitude. Our present 
model is geostrophically balanced on both the large and small scales, though as equation 
(12.501) shows a magnetostrophic balance does appear at higher, but subdominant, order 
in the mean momentum equation. A commonly employed dimensionless measure of the 
strong field dynamo is the Elsasser number, defined as the ratio of the Lorentz force to 
the Coriolos force and often expressed as 



2f Ip^r) 


(3.3) 


In terms of the nondimensional parameters and distinguished limits employed in the 
present study (e.g. Table Q]) we have 


A* = MPmRoRe = e. 


(3.4) 


Thus, as expected, the Elsasser number is small in our geostrophically balanced model. 
The above s caling is also consistent with previo us quasi-geostrophic magnetoconvection 
studies (c.f. I.Tones et~al. 2003; Gillet, et~aZl2007 ). In S74 and FC82 the Elsasser number 
is significantly smaller with A s = e 3 and A FC = e 2 , respectively. Although DNS studies 
are often characterized by order one Elsasser number, we note that these models are all 
limited to Pm ~ 0(1) and moderate values of the Ekman number. We speculate that if 
these models could increase their magnetic Reynolds number solely by increasing their 
Reynolds number and decreasing the Ekman number, rather than increasing Pm, they 
would begin to reach significantly smalle r Elsasser num bers; this trend may be apparent 
in the plane layer DNS investigation of Stellmach fe Hansenl (|2004r) and the spherical 
DNS study of Christensen fe Aubert ( 2006h . Also, we note that the value of the Elsasser 
number is also dependent upon how one scales the magnetic field. The form given by 
(EP1) is dependent upon the magnetic diffusivity, and independent of the velocity and 
length scales that are important for assessing dynamical balances; this is because it is 
typic ally assumed that the magnitude of the current density scales as J ~ UB/\xr\ (e.g. 
[Davidson 2001 1. To be consistent, the Elsasser number should be defined directly from 
the governing equations as 


_ —B ■ VB r2 

A _ pm _ ° (3 5) 

2fiz x u 2f Ip^lU 

In the present work we have 

A* = MRo = e 1/2 . (3.6) 

Similarly, for the Rm = 0(1) QGDM, this gives A* = e. This shows that although the 
Elsasser number is still small, it is independent of the magnetic diffusivity. Given that 
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the model developed in the present work is characterized by small Elsasser number, yet 
the magnetic energy is asymptotically larger than the kinetic energy for the low Pm 
QGDM, we can generally say that the partitioning of magnetic and kinetic energy is not 
necessarily indicative of dominant balances present within the governing equations. 

Finally, we note that it is possible to extend the CSDM to include fully nonlinear 
motions, a strong magnetic field, and small Pm by neglecting the large-scale horizontal 
modulations (i.e. d\ = dy = 0) appearing in equations (12.72M2.81I) : the result is given 
by 


dz^i/2 = ^ 0 Ol 


d T 0 o + dz(w' 0 6[) = —d 2 z 6 0 , 


<9 r B 0 = d z 


x (uo x ( b i/2>) 

^VlV’o + J(iPo, Wo) - d z< = z • V x (b 0 • Vx(b' 1/2 )) + Vitf 


Fm 


4 ^> 


dtw’o + Wo, O + Wo = ^Vi) + Bo • V 1 (&' 1 ( ; 2 ) ) + Vi w'o, 


Pr 


1 / 2 / 


d t (e[) + j(^, (9[)) + w' 0 d z e 0 = — v Koi), 

0 = B^-Viu' +^Lvi(b' 1/2 ), 

Fm 

9x{b'l0 + d y (b'$) = 0, 

and we note that equaton (12.761) is trivially satisfied in this case. 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


3.4. Applications to natural dynamos 

It is informative to compare the distinguished limits that were taken in the present work 
with what is known about the geodynamo and other natural dynamos; this exercise 
is useful for establishing the strengths and weaknesses of the present model. In the 
present work three timescales were necessary to allow time variations of the fluctuating 
convection, the mean magnetic field, and the mean temperature. The relative ordering 
of these timescales is given by 

t <C A t t -C A t T , (3.15) 

or 

t < e" 3/2 r < e~ 2 T, (3.16) 

upon noting that r = A~ x t and T = Ai^t. For instance, an order one change in the slow 
timescale r is asymptotically larger than an order one change in the fast timescale t. This 
ordering states that the large-scale magnetic evolution timescale lies midway between the 
small-scale convective timescale and the large-scale thermal evolution timescale. For the 
geodynamo, observations of the geomagnetic field show that t ~ O(10 2 ) years, r ~ 0( 10 4 ) 
years and T ~ O(10 9 ) years, suggesting that the above ordering is realistic. 

Utilizing the velocity based on the small-scale viscous diffusion timescale, the (small- 
scale) Rossby number can be related to the large-scale Ekman number as 

e = Ro = E]l 3 - (3-17) 

Studies suggest that the viscosity of the Earth’s core is similar to that of water at 
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standard temperature and pressure (e.g. Pozzo et al 120131 ) such that Eh ~ O(10 -15 ). 
Estimates for the small- and large-scale Rossby numbers, and the large-scale magnetic 
and hydrodynamic Reynolds numbers are given by, respectively, 


e = HT 5 , 

(3.18) 

Ro h = e 2 = 1CT 10 , 

(3.19) 

Rm H = RmE~ 1/3 = e~ 1/2 ~ O(10 2 ), 

(3.20) 

Re H = ReE~ 1/3 = 10 5 Re. 

(3.21) 


Although the small-scale Rossby number in the Earth’s outer core is not known, the above 
estimate for the la rge-scale Rossby number is not too different from estimates based on 
observations ( Finlav fe Amitl[2011 ). The large-scale magnetic Reynolds number possesses 
the correct order of magnitude estimate for the core. For the large-scale Reynolds number 
Ren we require knowledge of the currently unknown small-scale Reynolds number Re. 
Numerical simulations of the non-magnetic q uasi-geostro phi c co nvection equations can 
currently attain values of Re < O(10 3 ) (e.g. Ijulien et al\ 2012&h . leading to large-scale 
Reynolds numbers of Re h ~ O(10 8 ); this suggests the QGDM is in the appropriate 
dynamical regime necessary for understanding the geodynamo. 

The current model assumes an order one small-scale Reynolds number; coupled with 
the relationship Rm = RePm, this implies that Pm = O^ 1 / 2 ). The magnetic Prandtl 
number in the Earth’s core is thought to be Pm ~ O(10 -6 ), so it would seem that 
the distinguished limit taken for Pm is not quite in line with what typical values are 
in the Earth’s core. However, the theory developed here suggests that Pm = O^ 1 / 2 ) 
is sufficiently small to be in the asymptotic limit of small magnetic_Prandtl number; 
we can then reach smaller magnetic Prandtl numbers by reducing Pm. Furthermore, 
the distinguished limit of M = 0(e -1 / 2 ) shows that magnetic energy dominates when 
Pm <C 1; this result is thought to be generally consistent with what is known about the 
geo dynamo. 

As previously mentioned, we have chosen to non-dimensionalize the equations based 
on those dimensional quantities that possess the maximum amplitude; this implies that 
some non-dimensional parameters will contain a ‘mix’ of large- and small-scale dimen¬ 
sional quantities. Because of the importance of the small-scale motions, the dimensional 
length and velocity scales are then the small scales l and U , respectively. Pressure, tem¬ 
perature, and magnetic field were then based on the large-scale mean quantities. Any 
dimensionless parameter can be rescaled to be only a function of the small-scale or large- 
scale dimensional quantities by noting the asymptotic relations that were employed in 
the present work. For instance, the magnetic to kinetic energy ratio M is based on the 
large-scale magnetic field and the small-scale velocity field. When expressed solely in 
terms of the large-scale variables we have M = e _1 M « 10 7 where we’ve used the e 
estimate given by (13.1811 . 

Assuming a_large-scale magnetic field strength of O(10~ 3 ) Tesla for the core (e.g. 


g a l ar 

CJillet et aZ.ll201(Th . the traditionally defined Elsasser number (13.31) is A = 0(1). The model 
presented here is not capable of matching this value and represents one of the most signif¬ 
icant discrepancies between the present work and what is presently known about the geo- 
dynarno. With the exception of Jupiter, however, estimates suggest t hat the Elsasser n um- 
ber for other planetary dynamos is significantly less than unity ( Schubert & Soderlimd 
[201lh . As mentioned previously, the traditionally defined Elsasser number may not be 
an accurate indicator of the dominant force balance in the core as numerical simula¬ 
tions show that dynamos characterized by A = 0(1) remain geostrophically balanced 
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( Soderlund et al ■ 2012). With this result in mind we emphasize that some care must be 
taken with regard to estimating the dominant balance in DNS. Taking the curl of the 
momentum equation appears to be a common approach for determining the do minant 
force balance in simulations ( Christensen fc Aubert 20061 King fc Buffettl [2013h . Such 
a procedure can often lead to an erroneous conclusion, since the action of curling the 
momentum equation will automatically filter a dominant geostrophic balance. This is 
indeed the case for the present model where leading order geostrophic balances identi¬ 
fied in equations (12.4911 and (12.651) are removed by the curling operation. For instance, 
equation (12.521) shows that taking the curl of the mean momentum equation (12.491) yields 
the trivial result up to 0(e 1 ' 2 ). The first non-trivial result arises at 0(e) by curling the 
mean momentum equation (12.501) and will yield a vorticity equation involving a (vortic- 
ity) balance between the curl of the Coriolis and Lorentz force terms; this should not 
be interpreted as a force balance between the Coriolis and Lorentz forces. The salient 
point is that dominant balances can only be determined directly from the momentum 
equation, and not from the vorticity equation. Rather, analysis of the terms present in 
the vorticity equation for a geostrophically balanced flow will tell you which forces act 
as small perturbations to the main geostrophic balance. 

The alternative definition of the Elsasser number given by equation (|T5p . is O(10" 2 ) 
for the geodynamo, suggesting that the Lorentz force may not enter the leadin g order 
force balance for the large-scale dynamics of the core ( Nataf fc Schaeffeil 12015 ). More¬ 
over, previous work calculating the angular momentum of the Earth’s outer core must 
invoke geostrophy to obtain the large-scale velocity field within the core, and shows ex¬ 
cellent agree ment with values obtained independently from observations of the Earth’s 
rotat i on rate (Jault et al. 19881 : 1 Jackson et al. 19931 : Gillet et al. 20ld : Schaeffer fc Paisl 


2011! : Gillet et al. 2012 ). These results suggest that a large-scale thermal wind model, 


similar to the one developed here, may be relevant for understanding the physics of the 
geodynamo. Of course, the large-scale thermal wind model is characterized by Ax Az 
and is therefore unlike that of a deep spherical shell where the aspect ratio is order one. 
We were unable to obtain a closed large-scale model with Ax = that is coupled to 
the small-scale dynamics; future work is necessary to determine if this difficulty can be 
circumvented in a spherical geometry. 


The small-scale dynamics of the geodynamo are completely unconstrained due to lack 
of observations. The present work has utilized an expans ion base d on the smallness of the 
small-scale Rossby number. As discussed bv lNataf fc Schaeffer: (2015), the precise value 
of the small-scale Rossby number depends upon the detailed properties of the turbulence 
that is controlle d by t he re l ative importance of inertia, buoyancy and the small-scale 
magnetic field. iNataf fc Schaeffeil ( 20151 ) argue that although the Elsasser number, as 
given by equation (13.51) . is 0(0.01) for the largest scales of the Earth’s core, it may 
increase to 0(1 — 10) on the scale of about 10km such that the Lorentz force becomes 
more significant as the length scale is reduced. Indeed, as the QGDM shows, if the Lorentz 
force is subdominant on the largest scales then it must be present in the small-scale 
dynamics for the dynamo to reach a saturated, non-kinematic state. Future numerical 
simulations of the QGDM should therefore help to underst and t he s mall-scale dynamics 
of the geodynamo and other planetary magnetic fields (e.g. Stanley fc Glatzmaier 2010 ). 


We emphasize that exact agreement between the above estimates for the non-dimensional 
parameters and observations of natural dynamos is not required for the asymptotic the¬ 
ory presented here to accurately model convection-driven dynamos. As long as e is small, 
the QGDM represents an accurate reduced form of the full set of governing equations 
and can faithfully model geostrophically balanced dynamos in the plane layer geometry. 
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For atmospheric and oceanic flows, for instance, it is known tha t QG theory remains 
accurate for Rossby numbers just above e = 0(0.1) (e.g. Pedloskv 119871 ). 


4. Conclusion 

In the present work we have utilized a standard multiple scales asymptotic approach 
to develop a new multiscale dynamo model that is valid in the limit of small Rossby 
number. The small-scale model is characterized by a magnetically- modi fi ed ver s ion o f 
the quasi-geostrophic conve ction equations origi nally developed bv ljulien et ali 1 19981) 
and later studied in detail bv ISpraeue et al. ( 20061 ). The large-scale model is characterized 
by a thermal wind balance in the mean momentum equations; coupled with the mean 
heat equation, the large-scale model is equivalent to the well-known planetary geostrophy 
equations commonly employed for investigating the dynamics of the Earth’s atmosphere 
and ocean. We have discussed both low and order one magnetic Prandtl number models, 
showing that these two cases possess fundamentally different properties. For the low 
Pm case the magnetic energy dominates the kinetic energy, and ohmic dissipation is 
asymptotically dominant over viscous dissipation on the large-scales. When Pm = 0(1) 
it was demonstrated that the magnetic and kinetic energies become equipartitioned with 
weak large-scale ohmic dissipation. In both cases the dynamics are characterized by 
small Elsasser number since the motions are geostrophically balanced. The new model 
can be conside red a fully nonlinea r, gen eralized version of the dynamo model originally 
developed bv l Childress fc Sowardl ( 19721) . 

Numerical simulations of asymptotically reduced equation sets have proven useful for 
accessing dynamical regimes in rotating plane layer convection that are com put ation¬ 


ally demanding or i mpossible to reach with the use of DNS (e.g. ISorague et al\ 120061 : 


Julien et al. 2012 a\ H: Rubio et al. 2014 : Stellmach et aZ.l[2014h . These investigations have 
shown the dependence of the flow regime on both the Rayleigh and Prandtl numbers, 


new phenomena such as an asymptotic heat transfer sca ling regime (Julien et al. 2012qj). 


and large-scale vortex formation via an inverse cascade ( Julien et al 


20 1261 ). Future sim¬ 


ulations of the new model will help to shed light on the interaction of the se phenomena 
with a magnetic field. For instance, a recent DNS investigation bv lGuervillv et al] ( 20151) 
has shown that large-scale vortices can play a significant role in generating large-scale 
magnetic fields, particularly as the magnetic Prandtl number is reduced below one. 

Various extensions of the present work can also be carried out. While we have as¬ 
sumed, for simplicity, that the fluid is Boussinesq and that the gravity vector and rotation 
vector are aligned, it is a straightforward procedure to relax both of these constraints 
to include compre s sibility and varying ang le between the gravity and rotation vectors 
( Julien et a7ll200d : iMizerski fe Tobiasl l2013). Although the anelastic approximation ap¬ 
pears to agree well with the compressible equations for the case of order one Prandtl 


numbers ( Calkins et al. 20141) . it has recently been sho wn that it yields sp urious re- 


sults for low Prandtl number quasi-geostrophic convection ( Calkins et al. Il201~5h . However, 
Calkins et ali (|2015l) outli ned an app roach for extending the Boussinesq quasi-geostrophic 
convection equations of Sprague et al. I ( 20061) to the case of a fully compressible gas. 

The present model has focused on the plane layer geometry for the sake of mathemat¬ 
ical and physical simplicity. To further the applicability of the present work with that 
of natural systems, it will be useful to extend the present methodology to the rotat¬ 
ing cylindrical annulus ( Bussel[l986l) and to sp herical geometries. T he three-dimensional 
cylindrical annulus model recently developed bv ICalkins et ali ( 2013i) is particularly inter¬ 
esting since it possesses order one axial velocties, such that the Ekrnan pumping dynamo 
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effect investigated by iBusse ( 197 5) is no longer a prerequisite for dynamo action in this 


geometry. 


Acknowledgements 

The authors thank Stephen Childress, three anonymous referees, and Meredith Plumley 
for helpful comments that improved the manuscript greatly. This work was supported 
by the National Science Foundation under grants EAR #=1320991 (MAC, KJ and JMA) 
and EAR CSEDI #1067944 (KJ and JMA). 


REFERENCES 

Bender, C. M. & Orszag, S. A. 2010 Advanced mathematical methods for scientists and 
engineers I: Asymptotic methods and perturbation theory. New York: Springer. 

Busse, F. H. 1975 A model of the geodynamo. Geophys. J. Roy. Astr. Soc. 42, 437-459. 

Busse, F. H. 1986 Asymptotic theory of convection in a rotating, cylindrical annulus. J. Fluid 
Mech. 173, 545-556. 

Calkins, M. A., Julien, K. & Marti, P. 2013 Three-dimensional quasi-geostrophic convection 
in the rotating cylindrical annulus with steeply sloping endwalls. J. Fluid Mech. 732, 214- 
244. 

Calkins, M. A., Julien, K. & Marti, P. 2014 Onset of rotating and non-rotating convection 
in compressible and anelastic ideal gases. Geophys. Astrophys. Fluid Dyn. 000, 000. 

Calkins, M. A., Julien, K. & Marti, P. 2015 The breakdown of the anelastic approximation 
in rotating compressible convection: implications for astrophysical systems. Proc. Roy. Soc. 
A 471 (20140689). 

Chandrasekhar, S. 1961 Hydrodynamic and Hydromagnetic Stability. U.K.: Oxford University 
Press. 

Charney, J. G. 1948 On the scale of atmospheric motions. Geofys. Publ. 17, 3-17. 

Childress, S. & Soward, A. M. 1972 Convection-driven hydromagnetic dynamo. Phys. Rev. 
Lett. 29 (13), 837-839. 

Chini, G., Julien, K. & Knobloch, E. 2009 An asymptotically reduced model of turbulent 
Langmuir circulation. Geophys. Astrophys. Fluid Dyn. 103 (2-3), 179-197. 

Christensen, U. & Aubert, J. 2006 Scaling properties of convection-driven dynamos in ro¬ 
tating shells and applications to planetary magnetic fields. Geophys. J. Int. 166, 97- 114. 

Davidson, P. A. 2001 An Introduction to Magnetohydrodynamics. Cambridge: Cambridge Uni¬ 
versity Press. 

Dolaptchiev, S. I. & Klein, R. 2009 Planetary geostrophic equations for the atmosphere 
with evolution of barotropic flow. Dyn. Atmos. Oceans 46, 46-61. 

E, W., Engquist, B., Li, X., Ren, W. & Vanden-Eijnden, E. 2007 Heterogeneous multiscale 
methods: a review. Commun. Comput. Phys. 2 (3), 367-450. 

Fautrelle, Y. & Childress, S. 1982 Convective dynamos with intermediate and strong fields. 
Geophys. Astrophys. Fluid Dyn. 22 (3), 235-279. 

Favier, B. & Proctor, M. R. E. 2013 Kinematic dynamo action in square and hexagonal 
patterns. Phys. Rev. E 88 (053011). 

Finlay, C. C. & Amit, H. 2011 On flow magnitude and field-flow alignment at Earth’s core 
surface. Geophys. J. Int. 186, 175-192. 

Gillet, N., Brito, D., Jault, D. & Nataf, H.-C. 2007 Experimental and numerical studies 
of magnetoconvection in a rapidly rotating spherical shell. J. Fluid Mech. 580, 123-143. 

Gillet, N., Jault, D., Canet, E. & Fournier, A. 2010 Fast torsional waves and strong 
magnetic field within the Earth’s core. Nature 465, 74-77. 

Gillet, N., Schaeffer, N. & Jault, D. 2012 Rationale and geophysical evidence for quasi- 
geostrophic rapid dynamics within the Earth’s outer core. Phys. Earth Planet. Int. 202- 
203, 78-88. 

Grooms, I., Julien, K. & Fox-Kemper, B. 2011 On the interactions between planetary 
geostrophy and mesoscale eddies. Dyn. Atmos. Oceans 51, 109-136. 








22 


M. A. Calkins, K. Julien, S. M. Tobias and J. M. Aurnou 


Guervilly, C., Hughes, D. W. & Jones, C. A. 2015 Generation of magnetic fields by large- 
scale vortices in rotating convection. Phys. Rev. E 91 (041001). 

Haut, T. & Wingate, B. 2014 An asymptotic parallel-in-time method for highly oscillatory 
pdes. SIAM J. Sci. Comput. 36 (2), A693-A713. 

Jackson, A., Bloxham, J. & Gubbins, D. 1993 Time-dependent flow at the core surface 
and conservation of angular momentum in the coupled core-mantle system. Dynamics of 
Earth’s Deep Interior and Earth Rotation 72, 97-107. 

Jault, D., Gire, C. & Mouel, J. L. L. 1988 Westward drift, core motions and exchanges of 
angular momentum between core and mantle. Nature 333, 353-356. 

Jones, C. A. 2011 Planetary magnetic fields and fluid dynamos. Ann. Rev. Fluid Mech. 43, 
583-614. 

Jones, C. A., Mussa, A. I. & Worland, S. J. 2003 Magnetoconvection in a rapidly rotating 
sphere: the weak-field case. Proc. R. Soc. Lond. A 459, 773-797. 

Julien, K. & Knobloch, E. 2007 Reduced models for fluid flows with strong constraints. J. 
Math. Phys. 48 (065405). 

Julien, K., Knobloch, E., Milliff, R. & Werne, J. 2006 Generalized quasi-geostrophy for 
spatially anistropic rotationally constrained flows. J. Fluid Mech. 555, 233-274. 

Julien, K., Knobloch, E., Rubio, A. M. & Vasil, G. M. 2012a Heat transport in Low- 
Rossby-number Rayleigh-Benard Convection. Phys. Rev. Lett. 109 (254503). 

Julien, K., Knobloch, E. & Werne, J. 1998 A new class of equations for rotationally con¬ 
strained flows. Theoret. Comput. Fluid Dyn. 11, 251-261. 

Julien, K., Rubio, A. M., Grooms, I. & Knobloch, E. 20126 Statistical and physical bal¬ 
ances in low Rossby number Rayleigh-Benard convection. Geophys. Astrophys. Fluid Dyn. 
106 (4-5), 392-428. 

King, E. M. & Buffett, B. A. 2013 Flow speeds and length scales in geodynamo models: the 
role of viscosity. Earth Planet. Sci. Lett. 371, 156-162. 

Klein, R. 2010 Scale-dependent models for atmospheric flows. Annu. Rev. Fluid Mech. 42, 
249-274. 

Malecha, Z., Chini, G. & Julien, K. 2014 A multiscale algorithm for simulating spatially- 
extended Langmuir circulation dynamics. J. Comp. Phys. 271, 131 150. 

Miesch, M. S. 2005 Large-scale dynamics of the convection zone and tachocline. Living Reviews 
in Solar Physics 2 (1). 

Mizerski, K. A. & Tobias, S. M. 2013 Large-scale convective dynamos in a stratified rotating 
plane layer. Geophys. Astrophys. Fluid Dyn. 107 (1-2), 218-243. 

Moffatt, H. K. 1978 Magnetic Field Generation in Electrically Conducting Fluids. Cambridge: 
Cambridge University Press. 

Moffatt, H. K. 2008 Magnetostrophic turbulence and the geodynamo. In IUTAM Symposium 
on Computational Physics and New Perspectives in Turbulence, pp. 339-346. Springer. 

Nataf, H.-C. & Schaeffer, N. 2015 Turbulence in the core. In Treatise on Geophysics (ed. 
P. Olson & G. Schubert), , vol. 8. Elsevier. 

Parker, E. N. 1955 Hydromagnetic dynamo models. Astrophys. J. 122, 293-314. 

Pedlosky, J. 1987 Geophysical Fluid Dynamics. Springer-Verlag New York Inc. 

Phillips, N. A. 1963 Geostrophic motion. Rev. of Geophys. 1 (2), 123-176. 

Pozzo, M., Davies, C. J., Gubbins, D. & Alfe, D. 2013 Transport properties for liquid 
silicon-oxygen-iron mixtures at Earth’s core conditions. Phys. Rev. B 87, 014110. 

Roberts, P. H. 1988 Future of geodynamo theory. Geophys. Astrophys. Fluid Dyn. 44, 3-31. 

Robinson, A. & Stommel, H. 1959 The oceanic thermocline and the associated thermohaline 
circulation. Tellus 11, 295-308. 

Rotvig, J. & Jones, C. A. 2002 Rotating convection-driven dynamos at low Ekman number. 
Phys. Rev. E 66 (056308). 

Rubio, A. M., Julien, K., Knobloch, E. & Weiss, J. B. 2014 Upscale energy transfer in 
three-dimensional rapidly rotating turbulent convection. Phys. Rev. Lett. 112 (144501). 

Schaeffer, N. & Pais, M. A. 2011 On symmetry and anisotropy of Earth-core flows. Geophys. 
Res. Lett. 38 (L10309). 

Schubert, G. & Soderlund, K. 2011 Planetary magnetic fields: observations and models. 
Phys. Earth Planet. Int. 187, 92- 108. 



Multiscale dynamo model 


23 


Soderlund, K. M., King, E. M. & Aurnou, J. M. 2012 The influence of magnetic fields in 
planetary dynamo models. Earth Planet. Sci. Lett. 333-334, 9-20. 

Soward, A. M. 1974 A convection-drive dynamo: I. the weak field case. Phil. Trans. R. Soc. 
Lond. A 275, 611-646. 

Sprague, M., Julien, K., Knobloch, E. & Werne, J. 2006 Numerical simulation of an 
asymptotically reduced system for rotationally constrained convection. J. Fluid Mech. 551, 
141 174. 

Stanley, S. & Glatzmaier, G. A. 2010 Dynamo models for planets other than Earth. Space 
Sci. Rev. 152, 617-649. 

Steenbeck, M. & Krause, F. 1966 The generation of stellar and planetary magnetic fields by 
turbulent dynamo action. Z. Naturforsch. 21a, 1285-1296. 

Steenbeck, M., Krause, F. & Radler, K.-H. 1966 A calculation of the mean electromotive 
force in an electrically conducting fluid in turbulent motion, under the influence of coriolis 
forces. Z. Naturforsch. 21a, 369-376. 

Stellmach, S. & Hansen, U. 2004 Cartesian convection driven dynamos at low Ekman num¬ 
ber. Phys. Rev. E 70 (056312). 

Stellmach, S., Lisciiper, M., Julien, K., Vasil, G., Cheng, J. S., Ribeiro, A., King, 
E. M. & Aurnou, J. M. 2014 Approaching the asymptotic regime of rapidly rotating 
convection: boundary layers versus interior dynamics. Phys. Rev. Lett. 113 (254501). 

Stewartson, Iv. & Cheng, H. Iv. 1979 On the structure of inertial waves produced by an 
obstacle in a deep, rotating container. J. Fluid Mech. 91, 415-432. 

Welander, P. 1959 An advective model for the ocean thermocline. Tellus 11, 309-318. 



